Separation of Quasiparticle and Phononic Heat Currents in YBa 2 Cu 3 07 _$ 
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Measurements of the transverse (k xy ) and longitudinal 
(k xx ) thermal conductivity in high magnetic fields are used 
to separate the quasiparticle thermal conductivity (k xx ) of 
the CuC>2 — planes from the phononic thermal conductivity 
in YBa2Cu307_,5. k e xx is found to display a pronounced max- 
imum below T c . Our data analysis reveals distinct transport 
(t) and Hall (th) relaxation times below T c : Whereas r is 
strongly enhanced, th follows the same temperature depen- 
dence as above T c . 

PACS numbers: 74.72.Bk 



The study of heat transport in the superconducting 
state is well known to provide valueable information on 
the quasiparticle (QP) excitations and their dynamics. 
Compared to other probes of the QP-dynamics such as 
the microwave conductivity thermal transport has the 
advantage of probing only the QP-response, since the su- 
pcrfluid does not carry heat. On the other hand, a major 
complication in the analysis of the thermal conductivity 
is often a substantial phononic contribution k p x to the 
heat current. Such a situation is realized in the high-T c 
superconductors (HTSC), where theiieat current above 
T c is known to be dominated by fe^o. Accordingly, the 
interpretation of experimental data is ambiguous: For 
example, the maximum in the temperature dependence 
of the thermal conductivity in the superconducting state 
of YBa2Cri307_5 (YBCO) has been attributed to both a 
maxinuina of the electronic as well as the phononic contri- 
butionEl A clear separation of the phonon and QP heat 
currents is difficult and up to nowj-an unsolved problem. 

It has been pointed out in RcfsBEi that the transverse 
thermal conductivity k xy - the thermal analogue of the 
Hall effect - is free of phonons, i.e. k xy is purely elec- 
tronic. k xy is related to the electronic thermal Conduc- 
ted 
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tivity k xx according to k xy = k xx tan an where an is t 
thermal Hall angle. From the Wiedemann-Franz la 1 
one expects that an is equal to the electrical Hall angle 
an as obtained from the electrical (a xx ) and Hall (o~ xy ) 
conductivities via tan an — o~ xy /cr xx . In conventional 
metals tan an — oj c t where to c = eB/m is the cyclotron 
frequency and r is the usual transport relaxation time. 
In contrast, in the normal state of the HTSC tanevj/ 
is highly anomalous and has a temperature dependence 
which is distinctly different from t(T)H~II9. Up to now no 
data on the QP Hall angle are available below T c , but 
one may suspect that its behavior is anomalous as well. 
Therefore, in order to calculate k xx from k xy it is neces- 



sary to treat tan an = u c tr as an additional independent 
parameter which must be determined experimentally. 

In this letter we determine both k xx and tan an below 
T c in a single crystal of YBCO from combined measure- 
ments of k xy and k xx in high magnetic fields. Our main 
results are: (1) The electronic thermal conductivity of 
the Cu02-planes shows a pronounced maximum below 
T c which is strongly suppressed by a magnetic field. (2) 
tan an, as extracted from the thermal transport data be- 
low T c , displays the same temperature and magnetic field 
dependence as tan an obtained from electrical transport 
data above T c , and it passes smoothly through T c . This 
shows that tr ~ th and, remarkably, that th and r 
behave differently also below T c ! 

The thermal conductivity tensor k is the sum of an 
electronic and a phononic part k = k el + k ph . Jt-,is de- 
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fined via the heat current density = — /cVTliiJ. We 
assume that k ph remains diagonal even for B ^ 0, i.e. 
k P y = 0. In this case the transverse components of k are 
purely electronic. The transverse thermal conductivity, 
also called the Righi-Leduc-effcct, is measured as follows: 
In a magnetic field B = (0, 0, B) a temperature gradient 
V^T is applied in the a;-direction. Under the condition 
jh.y = 0a transverse temperature gradient V y T is found 
in the y-direction. Using the Onsager relations we find 
in this situation 
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with k-r 



k yy for twinned crystals without in-plane 



anisotropy. k xy can therefore be determined experimen- 
tally by measuring V x T, V y T, and the total longitudinal 
thermal conductivity k xx . 

Our measurements were carried out at constant tem- 
peratures with the magnetic field applied perpendicu- 
lar to the Cu02-planes. Typically temperature gradi- 
ents V X T of order 0.5 K/mm were applied using a small 
manganin heater mounted on top of the samples. The 
resulting transverse temperature gradients V y T of order 
1CU 3 K/mm in magnetic fields up to 14 T were mea- 
sured with AuFe-Chxnmel thermocouples calibrated in 
the same field rangeO. To eliminate offset voltages due 
to misalignment of the thermocouple we have measured 
for both field directions ±B in order to determine the 
Righi-Leduc component of V y T which must be antisym- 
metric with respect to field reversal. We have measured 
in 2 different modes: Either B was reversed at fixed tem- 
perature or we have heated the sample to temperatures 
above T c before the field was reversed. Because of vor- 
tex pinning effects this latter mode was used for all low 
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temperature measurements. We note that sweeping the 
magnetic field at fixed temperatures results in a strong 
hysteresis of k xx in almost the entire-temperature range 
below T c , similar as reported in RefH3; k xy , which must 
be extracted from the asymmetry of the field dependence, 
can therefore not be determined by sweeping the mag- 
netic field. We have tested our method by measurements 
on an insulator [k xy = 0) and on simple metalaU. Details 
of our experimental setup will be described elsewhere. 
The results presented here have been obtained on a high 
quality twinned single crystal of YBCO with dimensions 
1.9 mm x 2 mm x 0.38 mm and with a superconducting 
transition at T c ~ 90. 5K. 

Representative experimental results are shown in 
Fig. 1. k xx has a pronounced maximum at T max < T c 
which is strongly suppressed by the applied magnetic 
field. The absolute value of k xx (wlO W/Km at T c ), 
the relative upturn of k xx in zero field as characterized 
by the ratio k xx (T max ) / k xx (T c ) « 1.6 for our sample, as 
well as the sensitivity of the maximum, to magnetic fields 
are consistent with previous resultscl The overall tem- 
perature dependence of k xy is similar to that of k xx but 
the maximum of k xy occurs at higher temperatures and 
the relative change below T c is larger in comparison to 
k xx . The absolute magnitude af k xy is comparable to 
that reported previously in Ref.El. 
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Here k xx is the electronic contribution from the CuC>2- 
planes (Cu02-bilayers in YBCO) and k^, is that of the 
phonons. k xx describes a possible contribution from the 
CuO-chains which are present in YBCO along the b- 
direction of the orthorhombic structure. These chains 
are metallic for optimally doped samples and lead to 
a rather stroag a-6-anisotropy in untwinned crystals 
(crbb/cTaa — 2)c3. In a twinned crystal they should con- 
tribute to the electrical and the heat conduction on aver- 
age; a-6-anisotropy is of course absent. In the subsequent 
data analysis k c x h x must be treated differently from k xx 
since the CuO-chains as a quasi- 1-dimensional channel 
for charge and heat transport should neither contribute 
to the transverse transport coefficients a xy and k xy nor to 
the magnetic field dependence of the longitudinal ones. 
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FIG. 2. Ak xx = k xx (B) ~ k xx {B = 0) (•) and k xy /B (o) 
versus magnetic field B at fixed temperatures given in the 
figure. 

We compare in Fig. [2] Ak xx = k xx (B) — k xx (0) and 
k xy /B. Notably, these quantities have the same magnetic 
field dependence, i.e. 
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This observation provides the key to our data analysis. 
We define 



FIG. 1. Open symbols: k xx (upper panel) and k xy (lower 
panel) of YBa2Cu307_« as a function of temperature T for 
various fixed magnetic fields as indicated in the figure. Full 
symbols: k xx st as obtained from our data analysis (see text). 

For our data analysis we assume that in YBCO 3 chan- 
nels of heat conduction are present: 



k xy 

tan a R = lu c t r = j^j- 



(4) 



where tr(B,T) is a "relaxation time" introduced to 
parametrize the field and temperature dependence of 
tanafl. Using Eqs. @ and (Q) we find 
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Apparently, our experimental results suggest that the 
term in brackets is zero. Since the 3 terms in brack- 
ets refer to 3 distinct channels of heat conduction this 
requires that tr, k^ x , and k xx are separately field inde- 
pendent. Note that this is certainly reasonable in view 
of the origin of these contributions to the heat current. 

Assuming thus that the 3 terms in brackets vanish tr 
can be calculacted from our data by comparing k xx and 
k xy /B at different magnetic fields according to gtr / to = 
A{k xy / B) / Ak xx (see Eq. (g)). The result is shown in 
Fig. 3. Note that the values obtained for different mag- 
netic fields coincide within the experimental accuracy 
consistent with the anticipated field independence of tr. 

As a check of our result for tr we have also determined 
erff/m = o xy jBo xx for the same sample ptom measure- 
ments of a xy and a xx in the npcmal stateEJ. These data 
as well as their extrapolation!!^ to temperatures below 
T c are also shown in Fig. 3. The extrapolated values for 
t^ 1 look very similar to t^ 1 regarding the temperature 
dependence (see inset of Fig. 3), but they appear to be 
systematically larger by roughly a factor 2. However, 
note that tr as extracted from the thermal transport 
data is clearly unaffected by the CuO-chains and that 
this is also true for a xy . In contrast, a xx has a contribu- 
tion from the CuO-chains, i.e. a xx — o v xx + (cr xx ), where 
is the electrical conductivity of the Cu02-planes and 
(cr^) is an average of the chain contribution appropriate 
for a twinned crystal. With a"^ x w {<r xx p we conclude 
that eTjy/m = <J xy /Ba xx is underestimated by a factor 2. 
Correcting the normal state data for this factor we find 
excellent agreement between th and tr, i.e. our data 
tell tr ~ th- This strongly supports the procedure of 
our data analysis. 

Once Tfl is known the remainder of our analysis is 
straightforward: k e xx (B ^ 0) follows from Eq. ([|) using 
the data for k xy (B), and k xx st is obtained subsequently 
from Eq. (Q) for each field strength. As a test for inter- 
nal consistency we have verified that k xx st is indeed field 
independent. Finally, k xx (B = 0) follows from Eq. (|^) 
using k xx st and the zero field data for k xx . We have also 
determined k xx (B = 0) directly from Eq. (|4|) using an 
extrapolation of B/k xy to B — in good agreement with 
the results obtained from using k xx st and Eq. (||). 

In Fig. 4 we show k xx (B) as obtained from our data 
analysis. k xx represents the electronic thermal conduc- 
tivity of the Cu02-planes in YBCO. Our results thus 
confirm explicitly that k xx is strongly enhanced below 
T c . k xx oc TriQpT implies that r is strongly enhanced 
below T c overcompensating the decrease of the QP num- 
ber density uqp(T) with decreasing temperature. This 
confirms the results obtained for [the QP relaxation time 
from the microwave conductivity^. 




FIG. 3. T < T c : m/erji vs. temperature T obtained from 
eTR/m = A(k X y/B)/Akxx, 

where A(k xy /B) = k xy (Bi)/Bi - k xy (B2)/B 2 and 
Ak xx = k xx (B 1 ) — k xx {B2)- Different symbols correspond to 
different values of Bi and B 2 . T > T c : m/em obtained from 
a xy and a xx . Solid line: extrapolated normal state data (see 
text). Dashed line: extrapolated normal state data divided 
by a factor 2. Inset: The same data on a double logarith- 
mic scale. The dotted line corresponds to a T 2 temperature 
dependence. 

We also find that k xx is very sensitive to magnetic fields 
contrary to what is observed in the normal state, where 
the total thermal conductivity and thus k xx is field in- 
dependent. It is straightforward to attribute the field 
dependence below T c to an additional scattering mecha- 
nismO characteristic for the superconducting state such 
as scattering of QPs on vorticesDiiallj. Assuming the 
corresponding scattering rate r" 1 to be proportional to 
the number of vortices n v oc B one expects for the to- 
tal scattering rate that r" 1 = r^ 1 + t^ 1 = r^ 1 + &B. 
Here, Tj n includes the same scattering processes as in the 
normal state, i.e. it has in general an elastic defect and 
an inelastic contribution; the latter collapses below T c . 
Using k xy — k xx u c TR, = L(T)T<j xx u) c tr where L is the 
Lorentz-number and a xx = riQp(T)e 2 T /to we find 
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C(T)r- 1 - C(T) (r- 1 + aB) , 



(6) 



where C(T) 



2 1 \LTnQpe 2 TR) depends only on tem- 



perature. Thus a plot of B/k xy vs. B should yield a 
straight line. This is indeed the case as shown in Fig. 4. 
We note that by assuming a to be temperature inde- 
pendent Ti n (T) can be extracted from the slope and the 
intersection of the B/k xy vs. B curves (to within the 



constant factor a). 
below T r lii 



We find that Tj n increases strongly 
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FIG. 4. Left panel: Electronic thermal conductivity k x l x (B) 
of the Cu02-planes as a function of temperature for various 
magnetic fields as indicated in the figure. The zero field data 
have been obtained from k xx (B = 0) = k xx (B = 0)-k x % st (o) 
as well as from an extrapolation of B jk xy to B — (•). Right 
panel: B/k xy as a function of magnetic field B at various fixed 
temperatures given in the figure. 

kx% st as obtained from our data analysis is shown in 
Fig. 1. Remarkably, k x % st shows a pronounced maxi- 
mum below T c , too. Xhis maximum may be due to the 
phononic contribution!!!. However, it may also arise from 
the chain contribution. The latter has previously been 
determined experimentally from the a-6^aaisotropy of k xx 
in detwinned single crystals of YBCOcluJ. k xx shows a 
pronounced maximum below T c with an overall temper- 
ature dependence similar to that found here for fc£| st E9. 
Furthermore, the field independence of k^ x st implies that 
the phononic contribution k^ is independent of B. Such 
a conclusion has recently been drawn also on the basis of 
low temperature-results for the thermal conductivity in 
Bi-based HTSCES. 

Finally, our data show clearly that tr ~ th and that 
- as in the normal state - th and r behave differently 
also below T c . In particular, whereas r is strongly en- 
hanced below T c , th is unaffected by the superconducting 
transition and shows the same temperature dependence 
as above T c . This finding should provide important in- 
formation for the theoretical understanding of transport 
phenomena in the cuprates. 

In summary, we have presented a separation of the 
QP and phononic contributions to the thermal conduc- 
tivity below T c in YBCO based on measurements of the 
longitudinal and transverse thermal conductivity in high 
magnetic fields. Our data analysis shows explicitly that 
the QP contribution to k xx is strongly enhanced below 
T c and that it is the QP contribution to the heat current 
which is responsible for the magnetic field dependence of 
k xx ■ We find that - as in the normal state - 2 relaxation 
times must be distinguished also below T c : Whereas the 
QP relaxation time r is strongly enhanced below T c and 
magnetic field dependent, the Hall relaxation time th 
remains independent of B below T c and has the same 
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